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1. Introduction 

Operator splitting procedures are special finite difference methods one uses to 
solve partial differential equations numerically. They are certain time-discretization 
methods which simplify or even make the numerical treatment of differential equa- 
tions possible. 

The idea is the following. Usually, a certain physical phenomenon is the com- 
bined effect of several processes. The behaviour of a physical quantity is described 
by a partial differential equation in which the local time derivative depends on 
the sum of the sub-operators corresponding to the different processes. These sub- 
operators are usually of different nature. For each sub-problem corresponding to 
each sub-operator there might be a fast numerical method providing accurate so- 
lutions. For the sum of these sub-operators, however, we usually cannot find an 
adequate method. Hence, application of operator splitting procedures means that 
instead of the sum we treat the sub-operators separately. The solution of the orig- 
inal problem is then obtained from the numerical solutions of the sub-problems. 
For a more detailed introduction and further references, see the monographs by 
Hairer et al. [7], Farago and Havasi [6], Holdcn ct al. [8], or Hunsdorfer and Verwer 
[9]. 

Since operator splittings are time-discratization methods, the analysis of 
their convergence plays an important role. In our earlier investigations in Batkai, 
Csomos, Nickel [2] we achieved theoretical convergence analysis of problems when 
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operator splittings were applied together with some spatial approximation scheme. 
In the present paper we additionally treat temporal discretization methods of spe- 
cial form. Since rational approximations often occur in practice (consider e.g. Euler 
and Rvmgc-Kutta methods, or any linear multistep method), we will concentrate 
on them. Let us start by setting the abstract stage. 

Assumption 1.1. Assume that X is a Banach space, A and B are closed, densely 
defined linear operators generating the strongly continuous operator semigroups 
(T(t)) t >o and (S(t)) t >o, respectively. Further, we suppose that the closure A + B 
of A + B with D(A + B) D D(A) n D(B) is also the generator of a strongly 
continuous semigroup (U(t))t>o- 

For the terminology and notations about strongly continuous operator semi- 
groups see the monographs by Arendt et al. [1] or Engel and Nagel [5]. Then we 
consider the following abstract Cauchy problem 

[ ^- = (A + B)u(t), t>0, 

u(o) =xeX. 

For the different splitting procedures the exact split solution of problem (1) 
at time t > is given by 

< q (t) := [S(t/n)T{t/n)] n x (sequential), 

u^* (t) := [T(t/2n)S(t/n)T(t/2n)] n x (Strang), 

<(£) := [eS(t/n)T(t/n) + (I - Q)T{t/n)S{t/n)] n x with 9 € (0, 1) (weighted). 

However, in practice, we obtain the numerical solution of the problem (1) by 
$ applying a splitting procedure with operator A and B, 

^ defining a mesh on which the split problems should be discretized in space, 
and 

using a certain temporal approximation to solve these (semi-)discretized equa- 
tions. 

Thus, the properties of this complex numerical scheme should be investigated. In 
order to work in an abstract framework, we introduce the following spaces and 
operators, see Ito and Kappel [10, Chapter 4]. 

Assumption 1.2. Let X m , to G N be Banach spaces and take operators 

P m : X — > X m and J rn : X rn — > X 
having the following properties: 

(i) P m J rn = I m for all raefj, where I m is the identity operator in X m , 

(ii) lim J m P rn x — x for all x G X, 

m— >oo 

(hi) || J m \\ < K and ||P m || < K for all to G N and some given constant K > 0. 
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The operators P m , m £ N are usually some projections onto the spatial 
"mesh" X m , while the operators J m correspond to the interpolation method re- 
sulting the solution in the space X, but also Fouricr-Glcrkin methods fit in this 
framework. 

Let us recall the following definitions and results from [2]. First we assume 
that the exact solution u of problem (1) is obtained by using only a splitting 
procedure and discretization in space. 

Assumption 1.3. Let (A m , D(A m )) and (B m , D(B m )), m 6 N, be operators on X m 
and let (A,D(A)) and (B,D(B)) be operators on X, that satisfy the following: 

(i) Stability: 

Assume that there exist constants M > such that 

(a) \\Rc\R{\,A m )\\ < M and ||Re XR(X, A)\\ < M, 

(b) ||ReAi£(A, .B m )|| < M and ||Re XR(\, B)\\ < M 
for all Re A > 0. 

(ii) Consistency: Assume that P m D(A) C D(A m ), P m D{B) C D(B m ), and 

(a) lim J m A m P m x = Ax for all x G D(A), 

(b) lim J m B m P m x = Bx for all x G D(B). 

Til— >00 

Remark 1.4. a) If Assumption 1.3 is satisfied for M = 1, then by the Hillc 
Yosida Theorem A, B, A m , B m are all generators of contraction semigroups 
(T(t)) t > , {S{t)) t >o, (T m (t)) t > , (S m {t)) t > . Furthermore, from the Trotter^ 
Kato Approximation Theorem (see Ito and Kappel [11, Theorem 2.1]) it fol- 
lows that the approximating semigroups converge to the original semigroups 
locally uniformly, that is: 
Convergence: 

(a) lim J m T m (h)P m x = T(h)x Vx G X and uniformly for h G [0,£o]j 

(b) lim J m S m (h)P m x = S(h)x Vi £ X and uniformly for h G [0,£o] 
for any to > 0. 

(b) In turn, the resolvent estimates are satisfied if A, B, A mi B m are all generators 
of bounded semigroups, with the same bound for all m G N. One may even 
assume that these semigroups have the same exponential estimate, that is 

||r m (t)||, \\T{t)\\, \\S m (t)\\, \\S(t)\\ < A/e wt for all t > 0. 

In the following this would result in a simple rescaling that we want to spare 
for the sake of brevity. 

In order to prove the convergence of the splitting procedures in this case, we 
need to formulate a modified version of Chernoff's Theorem being valid also for 
the spatial discretizations. Our main technical tool will be the following theorem, 
whose proof can be carried out along the same lines as Theorem 3.12 in [2]. Let 
us agree on the following terminology. We say that for a sequence a m ^ n the limit 

lim a m , n =: a 

m,n— >oo 
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exists if for all e > there exists N G N such that for all n,m > N we have 
\\a m ,n - a\\ < e. 

Theorem 1.5 (Modified Chernoff-Theorem, [2, Theorem 3.12]). Consider a se- 
quence of functions F m : R + — > _Sf(X TO ), m G N, satisfying 

F m (0) = I m for all m G N, (2) 

and that there exist constants M > 1, w e R, such that 

\\[F m {t)] k \\x(x m) <Me k ^ for allt>0, m,keN. (3) 

Assume further that 

^ J™Frn{Jl)Prri£ JmPm% 

fc-vO ft 

uniformly in m G N, and £fta£ 

. . i . m -^m (ft) ^ Pm % / ~ \ 

Gx := hm hm (5) 

m— too h—>0 h 

exists for all x G D C A", where D and (Ao — G)Z? are dense subspaces in X 
for Ao > 0. Then the closure G of G generates a bounded strongly continuous 
semigroup (U(t)) t>Q) which is given by 



U{t)x= lim J m [F m m n P m x (6) 

m,n— >oo 

for all x G A uniformly for t in compact intervals. 



2. Rational approximations 

Our aim is to show the convergence of various splitting methods when combined 
with both spatial and temporal disretization. As temporal discretizations we con- 
sider finite difference methods, or more precisely rational approximations of the 
exponential function. Throughout this section, we suppose that r and q will be 
rational functions approximating the exponential function at least of order one, 
that is we suppose 

r(0) = r'(0) = 1 and g(0) = q'(0) = 1. 

Further, we suppose that these functions are bounded on the closed left half-plane 

C_ := {z G C : Rcz < 0}. 

Rational (e.g. A-stable) functions typically appearing in numerical analysis satisfy 
these conditions. An important consequence of the boundedness in the closed left 
half-plane is that the poles of r have strictly positive real part, and thus lie in 
some sector 

S e := {z : z G C, | arg(z) < 6} 
of opening half- angle 6 G [0, 
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It is clear that for an application of the Modified Chernoff-Theorem 1.5 
uniform convergence (w.r.t. m or h) plays a crucial role here (cf. [2]). Hence, the 
following lemma will be the main technical tool in our investigations. 

Lemma 2.1. Let A, A m , P mi J m be as in Assumptions 1.2 and 1.3. Let r be a 

rational approximation of the exponential being bounded on the closed left half- 
plane C_. Then for all x G D{A) we have 

J mX^J^A^y^P^X JmPmX 



h 

uniformly in m G N for h — > . 



JniAinr, Pjn % 



(7) 



The proof of this lemma is postponed to the end of this section. With its help, 
however, one can prove the next results: (1) on convergence of spatial-temporal 
discretization without splitting, (2) on convergence of the splitting procedures 
combined with spatial and temporal approximations. 

Theorem 2.2. Let A, A m , P m , J rn be as in Assumptions 1.2 and 1.3 and let A gen- 
erate the semigroup (T(t)) t >o- Suppose that r is a rational function approximating 
the exponential function bounded on C_ ; and that there exist constants M > 1 and 
w £ I with 

\\[r{hA m )] k \\ < Me kuh for allh>0, k,meN. (8) 

Then 

Urn J m r{±A m ) n P m x = T(t)x, 

m,n— >oo 

uniformly for t > in compact intervals. 

Proof. We will apply Theorem 1.5 with F m (h) = r(hA m ). The stability criteria 
(2)-(3) follow directly from r(0) = 1 and assumption (8). For the consistency (5) 
we have to show the existence of the limit in (4) uniformly in m G N. But this is 
exactly the statement of Lemma 2.1. □ 

Here is the announced theorem on the convergence of the sequential splitting 
with spatial and rational temporal discretization. 

Theorem 2.3. Let A, B, A m , B m , P m , J m be as in Assumption 1.1, Assumptions 
1.2 and 1.3, and let (U(t)) t >o denote the semigroup generated by the closure of 
A + B. Suppose that the following stability condition is satisfied: 

\\[q{hB m )r{hA m )] k \\ < Me khu} for all h>0, k, m G N. (9) 

Then the sequential splitting is convergent, i.e., 

lim [q(±B m )r(±A m )] n x = U(t)x. 

m,n— foo 

Proof. We apply Theorem 1.5 with the choice F m (t) = q(tB m )r{tA m ) for an arbi- 
trarily fixed t > 0. Since stability is assumed, we only have to check the consistency. 
To do that, first we have to show that 

hm = J m (A m + B m )P m x (10) 

/i->o h 
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for all x 6 D(A + B) and uniformly for radN. The left-hand side of (10) can be 
written as: 

Jrnq{hPrii^(fiA ra ^P ra X J m P m X 

h 

t / 1 T-t \ t-> Jm^iJ^-^-m^Pm-^ J m Pm , JrnQ{J^Pin)Pm^ ~~ JrnPm^ 

= J m q(hB m )P m 1 . 

h h 

Since the topology of pointwise convergence on a dense subset of X and the 
topology of uniform convergence on relatively compact subsets of X coincide 
on bounded subsets of S£{X~) (sec e.g., Engel and Nagel [5, Proposition A. 3]), 
it follows from Lemma 2.1 that the expression above converges uniformly to 

Jm(A m -\- B m ^ P m X . EH 

Theorem 2.4. Suppose that the conditions of Theorem 2.3 are satisfied, but replace 
the stability assumption with either 

\\[r{^A m )q(hB m )r('jA m )} k \\ < Me kh " for allh>0, k,meN 

for the Strang splitting, or 

|| [eq(hB m )r(hA m ) + (1 - e)r(hA m )q(hB m )f || < Me kh " 

for a O € [0, 1] and for all h > 0, k, m £ N in case of the weighted splitting. Then 
the Strang and weighted splittings, respectively, are convergent, i.e., 

lim [r(£A m )q(±B m )r{±A m )] n x = U(t)x (Strang), 
lim [Q q (±B m )r(±A m ) + (l-0)r{±A m )q(iB m )] n = U(t)x (weighted). 



n,ui— > oo 



Proof. The proof is very similar as it was in the case of the sequential splitting 
in Theorem 2.3. The only difference occurs in formula (10). In the case of Strang 
splitting we take 

F m (t) = ri^AM^BM^A^ 

and write 



h 

= J m r(±A m )q(hB m )P, ' 



+ t/ m T ( „ -A-m ) Pi 



t JmQ^iJ^PrrCjPm^ Jm Pm % t ^m^*( 2 -^-m ) Pm^ JmPm% 



By Lemma 2.1 this converges uniformly to 

Jrai^^^m ~t~ Pm ~t~ 2^ m ^FmX — Jm(A m -\- B rn )jP ra X. 

For the weighted splitting we choose 

F m (t) = eq{tB m )r(tA m ) + (1 - &)r(tA m )q(tB m ), 
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which results in 

J m [Qq(tB m )r(tA m ) + (1 - Q)r{tA m )q(tB m )]P m x - J m P m x 



e 



h 



+ (i-6) 



h 



h 

By using the argumentation for sequential splitting, this converges uniformly (in 
to) to 

QJm(A m + B m )P m x + (1 — @)J m (_B m + A m )P m x = J m {A m + B m )P m x 
as /i — !> 0. □ 

Proof of Lemma 2.1 

The proof consists of three steps, the first being the case of the simplest possible 
rational approximation, which describes the backward Euler scheme. The next two 
steps generalize to more complicated cases. 

Step 1. Consider first the rational function r(z) = — Then r(hA m ) = j^R(j^, A m ) 
for ft, > sufficiently small. Then the left hand side of (7) takes the form 

|| 7^ ( 7i < ^ m "^( 7i ' ^rn^jPmX J mPmX} J mAmPuiX 1 1 — 

— 1 1 7T ( 7i t ^ m "^( 7i ' ^ m ) ^tn ^ ( 7T -^-771 ) "^( X ' A 1n ) P m X^ JjnA^P^x]^ | — 

\\Jm (7^-^(7^ ^-m) ^m) A m P rn X || . 

For the uniformity of the convergence in to G N, take x G D(A). Then for A > 

XJ m R(X, A m ^ P m X — t/ m -R(A, Ajji^AmPmX -t- JmPmX . 

Hence. 

\\XJ m R(X,A m )P m x - J m P m x\\ < ||J mJ R(A,A m )P m || • \\J m A 

rri Pm X \ \ 

follows. By the stability in Assumption 1.3, 

K 2 A1 

\\J m R(X, A rn )P rn \\ < holds for A > 0. 

Further, by the consistency in Assumption 1.3 the sequence J m A rn P rn x is bounded. 
Therefore 

XJ m R(X,A m )P m x -)■ J m P m x (11) 

as A — > 00 uniformly in m G N. Since \\J m R(X, A m )P m \\ is uniformly bounded in 
to G N, we obtain by the densness of P m D(A) C X m that (11) holds even for 
arbitrary x G X. 

Since for x G D(A), by the consistency in Assumption 1.3, the set 
{ J m A m P m x : to G N} U {Ax} 
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is compact, for arbitrary e > there is N G N such that the balls B(JiAiPiX,e) 
for i = 1, . . . N cover this compact set. Now let m £ N arbitrary and pick i < N 
with \\JiAiPix — J m A m P m x\\ < e. Then we can write 

|| tt^tto-^Ct^) ^To) A^Pjx^x t/^^ylm-PTTixll 

— 1 1 ~i Jm-^i ji i A-yyi) A^PrnX i/i-A^-F^x -\- ^J m A m P m x «/j.A^.Fjx|| 

— II ~h'^ ril ^ J ^~h 5 ^'^O-^m {JmAmPmX J^A^PiX) || 
+ || {h^m-Ri^i A m )P m — I m ^jJ^Ai P{ X 1 1 + c 

f^Ce + || (j^J m R(j^, A m )P m — I m ) JiAiPix\\ + e, 

with an absolute constant C > being independent on to € N for h sufficiently 
small. Because of (11), the term in the middle also converges to as h — > 
(choosing A = 4). This proves the validity of (7) for our particular choice of the 
rational function r. 



Step 2. Next, let k G N and r(z) := (1 _ z 1 /fc)fc . 
r(hA m ) = [| A m )] fc . We have to prove 



Then r(0) = 1, r'(0) = 1 and 



-> (12) 



uniformly for m £ N as h — > 0. To do this we shall repeatedly use the following 
"trick": for y G D(A m ) we have y = R(t, A m )(¥ — -A m )y. Hence we obtain 

Jm-Pju-E — Jin T R ( Jl s -^-m ) J?nR{jii A m 'jA ni P m X 



Jtn\_hRi.h^ A-Tn ) ] *^ 



fe-1 
3-0 



By inserting this into the left hand side of (12) we get 



J-m ( Jl R(, "fa : ) ) ^rr 



(13) 



1 

3=1 

By Step 1, we have J m 4i?(-r, A m )]A m P m x —> J m A m P m x uniformly in m 6 N as 
/i — > 0. Now, since J m [|i?(| , Am)} 3 P m is uniformly bounded for m G N, 1 < j < fc 
and h small enough, we also have that J m [j-R(j-, A m )] 3 A m P m x — > J m A m P m x 
uniformly in to G N as h — > for each 7 = 1, . . . k. This shows that the expression 
in (13) converges to uniformly in to G N. 



Step 3. To finish the proof for the case of a general rational function 



Splitting with approximations 



9 



we use the partial fraction decomposition, i.e., we write 



•c*) = EEti 



C 



with some uniquely determined Cy € C. Since, by assumption, r(0) = 1 and 
r'(0) = 1, we obtain 

i=l j=l i— 1 j — 1 

Since r is bounded on the left half-plane, we have that the poles A, of r have 
positive real part, Re Aj > 0. For j = 1, . . . , i/j, i = 1, . . . , / consider the rational 
functions 

nj(z) := 



1 



(i-*/i) J " 

and the operators Aij <m := j-A m . Then 

r ( z ) = EE C ^ T *» 

i=l 3=1 

We shall apply Step 2 to these rational functions and to these operators. To do 
that we have to check if the required assumptions are satisfied. Obviously, have 
the properties needed. The consistency part of Assumption 1.3 is trivially satisfied 
for Aij tm , P m and J m , m S N. The uniform boundedness of XR(X,Aij iTn ) follows 
from the stability Assumption 1.3. (a). Indeed, that condition implies 

||AiZ(A, A m \\ < M<f, for all A e S^, 

where is any sector with opening half-angle <f> € [0,5-). If Aj £ £ v is a pole of 
r then £ Eg for A > 0. So have 



\XR(X,±A 



\Xi r>l\\i 



\\XR{X,Aij, m )\ 
By Step 2, we have that 



} R(±f,A m )\\ <M 9 for all A > 0. 



JmAij^ffiPjyiX y 



uniformly in m € N as h — > 0. By taking also the equalities (14) into account this 
yields 

1 r 



Jm7'(hAj n ^P rr iX JrnRmX 



J yji A rn Prn ^ 



i=i i=i »=i i=i 

^ ^ ^ ^ T~ (^rn^ij{h~^A})P rtlj X — J m P 7n X^ — Jm A-ij _ m P 7n 



i=i j=i 
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uniformly in m £ N as h — > 0. This finishes the proof. □ 

Finally we remark that in the present paper we only treated an autonomous 
evolution equation (1). In the case of time-dependent operators A(t) and B{t) 
we have already shown the convergence in [3] for numerical methods applying 
splitting and spatial discretization together. The extension of our present results 
concerning the application of an approximation in time as well, will be the subject 
of forthcoming work. 
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